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1. Introduction 

In the last few years lattice calculations have stalled achieving the level of precision and the 
control of uncertainties necessary to extract phenomenologically relevant results |[T]]. Simulations 
with dynamical quarks are required for this task, since the uncontrolled errors associated with 
the quenched approximation are usually the main source of uncertainty in these calculations. An 
example of this fact can be seen in the study of indirect CP violation in the neutral kaon system. 

The CP violating effects in K° - K° mixing are parametrized by Ek- Experimentally, this 
quantity is known with a few per-cent level precision. On the other hand, theoretically it is given 
by the hadronic matrix element between K° and K° of the AS = 2 effective hamiltonian 

H$f = C A 5= 2 (aO / d 4 xQAs= 2 (x) (1-1) 

with 

Gas=2 (x) = [saYfid a ] yA (x) [spfdp] V _ A (*)• (1-2) 

The Wilson coefficient Cas=2(a0 * s a perturbative quantity known to NLO in a s in both NDR 
and HV schemes. It depends on several CKM matrix elements about which we would like to obtain 
information. The matrix element (K°\Qas=2\K°), that encodes the non-perturbative physics of the 
problem, is usually normalized by its VIA value, defining Bk as the ratio 

B (u) (K°\Q AS =2^)\K°) a 3) 

KW l(K°\s Y ,Y 5 d\0)(0\sY,Y 5 d\K°)- 

The main source of uncertainty when one tries to constrain the value of the combination of 
CKM matrix elements involved in the theoretical calculation of Sk using its experimental value is 
the error associated to the determination of B K [ph. Improvement in the calculation of B K is thus 
crucial in order to get information about the unitarity triangle. 

The value of Bk that the phenomenologists are using at present in their studies of the unitarity 
triangle that is considered as the benchmark of the lattice calculations of this parameter, was 
obtained by the JLQCD collaboration [||] using unimproved staggered quarks in the quenched 
approximation. The value given in H is B N K DR {2 GeV) = 0.628(42). 

The main source of uncertainty in this calculation is the unknown error from quenching, that 
could be as large as a 15% according to the ChPT estimate performed in [Q]. In order to have a 
prediction at a few per-cent level it is thus necessary to perform a dynamical calculation of Bk 
that eliminates the quenched uncertainties. Another drawback of the calculation in ^ is the fact 
that it is affected by large scaling uncertainties. We will see that the scaling behaviour is going to 
be much better using improved 1 staggered actions instead of the standard unimproved staggered 
action used by the JLQCD collaboration. A third way of improving the JLQCD calculation would 
be incoiporating SU (3) breaking effects by using kaons made up of non-degenerate quarks, instead 
of degenerate quarks with m s /2. 

The goal of this work is to perform a dynamical calculation of Bk that eliminates the quenched 
uncertainty, using improved staggered fermions that have been proved to reduce the large iff (a 2 ) 
discretization errors generated by the taste-changing interactions. 

1 In the improved actions the thin links are substituted by fattened links and the quark-gluon interactions that violate 
the taste symmetry are reduced. 
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p 


Volume 


Wconfs 


fl _1 (GeV) 


m s /2 


5.7 


12 3 x 24 


150 


0.837(6) 


0.086/0.064 


5.93 


16 3 x 32 


50 


1.59(3) 


0.039/0.030 



Table 1: Parameters in the quenched simulations. The values of m s /2 are for the HYP and Asqtad staggered 
actions respectively. 
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Figure 1: Scaling of Bj[ (2GeV) with a for improved staggered actions compared to the JLQCD unim- 
proved staggered results. 

Preliminary results for that study were presented in The matching coefficients needed in 
the calculation of the renormalized with the action used in our dynamical simulations were not 
available at that moment, so an approximate renormalization was performed in order to get those 
preliminary results. The correct renormalization coefficients were calculated later [|]] and have 
been used to obtain the results reported in these proceedings. 

2. Improved versus unimproved staggered actions: scaling violations 

The first thing we analyse is the impact of using improved staggered actions in comparison 
with the unimproved staggered action used in [||], that suffers from large scaling violations. We do 
this analysis for two different improved actions: the HYP [Q] and the Asqtad 2 

The values of the parameters used in the simulations are shown in Table [T]. We choose these 
parameters to be the same as those used by the JLQCD collaboration in order to make a clear 
comparison with their results. In particular, we match kaon masses at a given j8 to those of the 
JLQCD collaboration. The matching of the lattice operators to the continuum ones in the MS 
scheme have been performed perturbatively using the one-loop coefficients calculated in [Q]. 

The results we obtain for 5^ DR (2GeV) as a function of the lattice spacing are shown in Figure 
|TJ In that Figure a clear improvement in the scaling can be seen when using improved actions, 

2 In our quenched analysis we are using a variation of the Asqtad action with no improvement in the gauge action. 
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p 


ftcanfs 




m s /2 


a v (l/a) 


6.76 


560 


0.01/0.05 


0.02 


0.4723 


6.79 


414 


0.02/0.05 


0.02 


0.4699 



Table 2: Parameters in the dynamical simulations. 



in particular in the HYP case. We expect such improved scaling to survive unquenching and this 
therefore allows us to perform reliable dynamical calculations with only a few values of the lattice 
spacing and even obtain valuable information from a single point simulation. 



3. Dynamical value of Bk 

We now incorporate dynamical effects in the calculation of Bk using one of the improved 
staggered actions analyzed in the quenched approximation, Asqtad, since the final goal is to elim- 
inate the irreducible systematic error associated to quenching that dominates the final uncertainty 
in previous determinations of Bk- 

We performed a dynamical calculation of Bk with the Asqtad action, using the configurations 
from the MILC collaboration with nf = 2+1 dynamical flavours [Q]. The results reported in 
these proceedings correspond to the analysis at one lattice spacing with a = 0.125 fm and two 
different values of the sea light quark masses. The parameters used in the dynamical simulations 
are collected in Table ||| 

The conversion of the values of the bare lattice operators to a value for B ! ^ DR (2Ge\) has 
been done perturbatively using the &(cc s ) lattice to continuum matching coefficients in [Qj. In the 
matching process we take a s in the V scheme at the scale 1 ja where a is the lattice spacing -see 
values in Table ^[ The value for GCy with Nf = 3 has been taken from the recent 4-loops lattice 
determination in [|To|]. 

The results we obtain for B K J)R (2GeY) including only statistical errors as a function of the 
light sea quark masses are shownin Figure |2[ A decrease of the value of Bk with the reduction 
of the dynamical quark mass can be appreciated in that Figure. After linearly extrapolating these 
results to the physical s and u(d) masses following the staggered chiral perturbation study in [|Tl|], 
the value of B N K DR (2GeV) we obtain is 

B K DR (2GeV) =0.618(18)(19)(30)(130), (3.1) 

where the first error is statistical, the second is from chiral fits, the third one is from discretizations 



errors and the final one is from the perturbative conversion to the MS scheme. The value in ( p.l| ) 
is equivalent to B K = 0.83 ±0.18, with B K defined as the product of S^ DR (2GeV) and the Wilson 
coefficient in the effective hamiltonian ( |lTl| ) in the MS — NDR scheme and at 2GeV. This quantity 
is scheme and scale independent at ff(ot s ). Note that this value of Bk is very similar to the previous 
quenched staggered result in [^] (Bk = 0.86(6)(14)), so any final conclusion about the enhance- 
ment or decreasing due to the inclusion of dynamical effect in its calculation would need a drastic 



reduction of the error quoted in (3.1). That error is dominated by the uncertainty associated to the 



possible &(ccf) corrections in the lattice to continuum matching process, that we estimate to be of 
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Gauge invariant B (2GeV): dynamical vs. quenched 
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Figure 2: Dynamical value of B^ DR (2GeV) as a function of the ratio between the light sea quark mass and 
the (real) strange quark mass. The lines represent the quenched results from [Q]. 

the order of (ay(l/a)) 2 . At the lattice spacing we are working, the unquenched values of ay are 
given in Table ^ and translate into a ~ 20% error in the result for Bk- 



4. Summary and discussion 

Most of the previous lattice calculations of Bk, in particular that used in the unitarity triangle 
analysis, were performed in the quenched approximation. This induces a large, essentially un- 
known and irreducible systematic error into the result. Precise simulations with dynamical fermions 
are necessary in order to be able to make full use of the experimental data on Ek to constrain the 
CKM matrix. These dynamical simulations are feasible with present computers using staggered 
fermions at light dynamical quark masses. However, the unimproved staggered action suffers from 
large taste-changing interactions that generate important scaling corrections, as those found by the 
JLQCD results. We have shown in Section [2| that the scaling behaviour is much better when us- 
ing improved staggered actions, that have been designed to reduce the non-physical taste-changing 
interactions as well as other cut-off effects. 

As a first step in the dynamical study of Bk on the MILC configurations we have calculated 
this quantity in two ensembles at a = 0. 125 fm and with two different light sea quark masses. 
In doing that, we have used the recent results for the matching coefficients corresponding to the 
Asqtad action [0], the one used in the numerical simulations. The value of the renormalized Bk 
in the MS scheme we obtain is B% DR (2GeV) = 0.618 ±0.136 or, equivalently, B K = 0.83 ±0.18. 
No sizeable deviation from quenched results can be inferred from this preliminary value. The 
uncertainty in this determination is dominated by the perturbative error associated to the one-loop 
matching process, so, in order to have the precision of a few per-cent needed by phenomenology, a 
two-loop matching or a numerical matching method is required. 



The next step will be to redo our calculation on a finer lattice [12]. This new result for a 
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different lattice spacing will allow us to perform an appropriate continuum extrapolation 3 , reducing 
the final discretization errors. The perturbative error will also be reduced since a,y(\/a) is smaller 
for finer lattices. And we will be able to check whether the good scaling behaviour observed in 
the quenched approximation using improved staggered actions (Asqtad) is in fact present in the 
dynamical results. 

We also want to investigate other issues in the study of Bk, like the differences between using 
invariant and non-invariant gauge operators, the application of the staggered chiral perturbation 
theory results in [O], the impact of SU (3) breaking effects or the chiral limit value of this quantity 
that could be compared to recent continuum calculations [p"3|]. 
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In view of the improved scaling behaviour we obtained in the quenched approximation within the Asqtad action 
described in Section 0, results for two lattice spacings will be enough to perform a reliable continuum limit. 
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